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Recent experiments show oscillations of period h/2e in conductance vs. magnetic flux of charge 
density wave (CDW) rings above 77 K, thus revealing macroscopically observable quantum behavior. 
The time-correlated soliton tunneling model discussed here is based on coherent, Josephson-like 
tunneling of microscopic quantum solitons of charge 2e. The model interprets the CDW threshold 
electric field as a Coulomb blockade threshold for soliton pair creation and draws upon the theory 
of time-correlated single-electron tunneling to interpret CDW dynamics above threshold. As in 
Feynman's derivation of the Josephson effect, the picture treats the Schrodinger equation as an 
emergent 'classical' equation to describe the time-evolution of Josephson-coupled fluidic complex 
order parameters. Coupling of vector or time- varying scalar potentials to the order parameter 
phases enables h/2e quantum interference in CDW rings and also leads to interesting behavior of 
CDWs in response to oscillatory electric fields. The ability to vary both magnitudes and phases of 
the Josephson-coupled order parameters is crucial to any future applications in quantum computing. 



I. INTRODUCTION 

Recent developments necessitate a transformation in 
our understanding of charge density wave (CDW) trans- 
port to one based on quantum principles P The CDW is 
a correlated electron (or electron-phonon) system that, 
like a superconductor, can transport electrons through a 
quasi-one-dimensional or layered crystal en massed It is 
the only known such system capable, in the linear chain 
compound NbS3, of collectively carrying electric current 
above 37° C, the temperature of the human bodyP More- 
over, a significant body of evidence highlights the im- 
portance of CDWp^ stripe and other charge- and/or 
spin-ordered phases in high-T c and other unconventional 
superconductors as carrier concentration is varied, e.g. 
by doping. Some experiment^ suggest possible interfa- 
cial superconductivity or a related phase transition near 
the boundary between ion-implanted and unimplanted 
regions of a CDW in NbSe3. 

Importantly, Aharonov-Bohm (A-B) quantum interfer- 
ence effects, showing oscillations of period h/2e in CDW 
conductance vs. magnetic flux, have been reported^ for 
TaS3 rings above 77 K. Similar h/2e oscillations have 
been reproduced, as reported in 2012, 11 for at least five 
TaS3 rings with circumferences of up to 85 jam. The A-B 
oscillations are only observed in the CDW (not normal 
electron) magneto-conductance above the threshold elec- 
tric field for CDW transport, and the amplitude of the 
oscillations scales with CDW current. The ring experi- 
ments show that, at least for some materials, the CDW 
condensate exhibits quantum phase coherence over dis- 
tances of several ten's of microns. There is no classical 



sliding CDW explanation for such extraordinary behav- 
ior, which manifests Planck's constant at the macroscopic 
scale. These results thus underscore the need for a fun- 
damental paradigm shift, in which the laws of quantum 
physics play a crucial role in describing CDW electron 
transport. 

One ring was reported^ to exhibit telegraph-like tem- 
poral switching between high and low CDW current 
states, the high current state showing substantially larger 
amplitude h/2e oscillations than the low current state. 
This telegraph-like near destruction and reappearance 
of A-B quantum interference in the CDW ring indicates 
quasi-periodic loss of quantum coherence, suggesting two 
types of transport involving either probabilities or prob- 
ability amplitudes (the former lacking vs. the latter in- 
cluding quantum coherence). Alternatively, the behav- 
ior may suggest some form of macroscopically observable 
"wave- function collapse." Regardless of which interpre- 
tation emerges, a deeper understanding of the observed 
behavior based on the laws of quantum physics could ul- 
timately prove important to condensed matter physics 
and possibly to the foundations of quantum physics. 

Earlier experiments have shown that, in NbSe3 and 
TaS3, the CDW displaces very little below the thresh- 
old electric field for CDW transport. This is evident in 
NMR experiments^ showing a 2° CDW phase displace- 
ment in NbSe3, as compared to the classically predicted 
90° displacement just below threshold. Further evidence 
is provided by dielectric and other ac response (mixing, 
etcP^Qj) measurements, which exhibit a flat bias depen- 
dence as compared to the classically predicted divergent 
dielectric response shown in Fig. [TV a). These experiments 
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reveal that, even just below threshold, each portion of 
the CDW sits near the bottom of a pinning potential 
well. This suggests that, for these samples, the mea- 
sured threshold is substantially smaller than the classical 
depinning field and likely a Coulomb blockade threshold 
for charge soliton nucleation.^I3E3 
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FIG. 1. (a) Bias dependent dielectric response, showing classi- 
cal predictions vs. experiment. Classical models include clas- 
sical sine-Gordon, incommensurate harmonic chain (CF), 19 
random pinning,^ renormalization group (NM) 21 models, 
and where f = 1 — E/Et> Some NbSe3 measurements were 
carried out in our lab using a bridge circuit, while additional 
measurements were carried out by zcP and WMG 23 at the 
indicated temperatures and frequencies, (b) Experimental 
CDW conductance vs. electric field for NbSe3 as compared to 
the Zener tunneling curve exp[— Eo/E] (solid line) pointed out 
by BardeenP^ Adapted with permission from ref. ^H, Copy- 
right 1990, American Institute of Physics. 



The I-V characteristic has been found, starting with 
the early experiments on NbSe3p^ to progress from a 
rounded Zener tunneling-like char act eristic^ (Fig. [ljb)) 
to a nearly piecewise linear form in crystals with fewer 
impurities. This behavior further supports soliton pair 
creation with a Coulomb blockade threshold! 1 * 17 * 181 Soli- 
ton pair creation is analogous to Landau- Zener tun- 



neling, recently applied to Fe rmi superfluid gases 
Schwinger pair production] 27 1 28 1 or creation of supercon- 
ducting vortex-antivortex pairsP^ The existing classical 
models fail to explain the shape of the CDW /- V curves 
of NbSe3 and TaS3 m a straightforward fashionPS Fi- 
nally, there is no existing classical sliding interpretation 
for t he h/ 2e quantum interference effects seen in CDW 
rings ESE] Any viable CDW transport theory of this ex- 
traordinary phenomenon must contain Planck's constant, 
even at the macroscopic level. However, this does not 
rule out the possibility of using the Schrodinger equation 
itself as an emergent 'classical' equation, as discussed by 
Feynman in the context of superconductivity (refP^, vol. 
Ill, Ch. 21). This approach, novel for CDWsP of em- 
ploying the Schrodinger equation to describe classically 
robust complex order parameters will be discussed in sec- 



tion III The following section discusses a modified sine- 
Gordon model, the simplest possible model of a pinned 
CDW. 



II. PINNED CHARGE DENSITY WAVE AS 
MASSIVE SCHWINGER MODEL 

A CDW has a modulated charge p(x,t) = po(x,t) + 
pi cos[2kFX — c/)(x, t)] along the axis of a linear chain com- 
pound. Here po(x,t) contains background charge and 
any excess or deficiency of charge oc d(f)/dx. The entire 
CDW condensate and Peierls gaps, initially at ±fcp, can 
be displaced in momentum space, resulting in a current: 
hdw oc d(j)/dt^^ Although a real CDW is pinned by 
impurities, in some materials it will still transport a cur- 
rent provided the applied field E exceeds a threshold Et- 
Displacing the CDW by one wavelength (advancing </> by 
2tt ) returns the system to its original state (except for 
charge displaced between contacts) so the pinning energy 
is periodic in <fi: u p [l — cos^]. (A qu antum version^ of 
the Fukuyama-Lee-Rice (FLR) model ESIISD including dis- 
order, is more accurate but observed voltage oscillations 
suggest the simple sine-Gordon picture captures much of 
the physics for high quality crystals.) 

Unlike a superconductor, the CDW charge modulation 
does not couple directly to a uniform electric field or vec- 
tor potential. However, gradients or kinks in phase carry 
charges that (1) couple to an externally applied field and 
(2) generate their own electric fields that lead to electro- 
static interactions. These electrostatic interactions be- 
tween kinks, often neglected in previous theories, are im- 
portant whether treating the system classically or quan- 
tum mechanically. If the CDW phase is initially fixed at 
zero at the contacts or at ±oo, advancing the phase by 
(j) in the middle creates charged kinks that produce an 
internal field: E^ = (E* /2tt) , where E* = 2e/(eA ch ) 
is the field created by a 2tt phase soliton-antisoliton pair 
and A c h is the cross-sectional area per chain. 

Figure [2] shows the combined effects of the applied field 
E and the field E* created by a pair of soliton domain 
walls. The difference in electrostatic energy densities, 
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\e(E ± E*) 2 — \eE 2 , with and without the pair is posi- 
tive when E is less than the Coulomb blockade threshold 
field, Et — \E* — en c h/€. Here n c h = 1/A c h is the num- 
ber of parallel chains per unit cross-sectional area. The 
empiricall y obs erved relation eE T ~ en c h pointed out 
by Grimer^ 3 -^ thus emerges naturally from this picture. 
The simplest classical model predicts^ 3 -^ the classical 
depinning field E c \ to scale as: eE c i = 47ren c h, where 
E c i oc n 2 for weak pinning (m being the impurity con- 
centration). This yields: Et = E c i/4tt, which has the 
same impurity dependence as E c \ for a fixed tempera- 
ture. Screening by normal carriers further enhances e 
and reduces the ratio: E T /E c i. For fixed n^, the tem- 
perature dependence of carrier concentration and e leads 
(inversely) to the strong temperature dependence of Et 
seen in some materials. In addition, en c h is multiplied by 
the condensate fraction p c in a more precise description. 



(a) 




FIG. 2. (a) CDW phase vs. position, showing internal field E* 
produced by a soliton-antisoliton domain wall pair, (b) Model 
of density wave capacitance showing nucleated domain walls, 
more realistically depicted as soliton dislocation droplets in 
section [Hi] The applied field E partially or completely cancels 
the internal field E* . 

Following the quantum field theory literature, E re- 
lates to the 'vacuum angle' as 6 = 2tt(E/E*). For phase 
displacements <\> between contacts, E partially cancels 
E^, yielding an electrostatic energy ue{0 — (f)) 2 J^jThe po- 
tential energy per chain can then be written a d 1 * 18 * 39 !: 

U[<j>} = j dx [2u p [l- cos </)(x)]+ue(0 - <K^)) 2 } (!) 

This is a variant of the bosonic massive Schwinger model, 
studied as a model of quark confinement! 1 6 * 4Q I and first 
adapted to explain the CDW quantum threshold field by 
Krive and RozhavskyP^ The usual linear coupling oc — Ocj) 
is contained in the quadratic term, as are electrostatic 
contributions oc (j) 2 and 2 . When < the system 
is stable classically and quantum mechanically (Fig. [3| . 
When > 7r, the (j) ~ 2ir state becomes the lowest energy 
state. 




FIG. 3. (Left) Potential energy vs. for two values of 
with many degrees of freedom illustrated as blue dots. Tun- 
neling can only occur if > ty(E > Et), when "bubbles" of 
the phases (j)k for the parallel CDW chains can nucleate by 
tunneling into the adjacent well. (Right) Potential energy 
parabolas u vs. 0, in which the phases (j)k ~ 4> are sitting 
in various potential minima, (j) ~ 27m. The first crossover 
between parabolic branches occurs at = tt. 

Thus, = <j) demarcates the boundary 39 above the 
system can decay into the lower well. Several quasi- 1-D 
systems appear to be in the sweet spot of interchain inter- 
actions - strong enough to avoid being swamped by ther- 
mally excited soliton dislocations but not strong enough 
to remain forever trapped in the higher well. Some NbSe3 
crystals suddenly switch into a higher CDW current- 
carrying state as the field is increased^ and show a hys- 
teretic I-V curve. A natural interpretation is that, as 
is increased above 7r, the system is temporarily trapped 
in the higher metastable well (Fig. [3]) before decaying 
rapidly into the lower well. Other materials show more 
than one threshold field SHHl The picture here provides a 
simple interpretation: that the lower threshold field is the 
Coulomb blockade threshold for soliton nucleatior i 1 * 17 * 18 * 
while the upper threshold is the classical depinning field. 

Figure [ZFa) shows blue bronze datsP^l that, especially 
at 48 K, exhibits two distinct threshold fields above which 
the conductance increases. The upper threshold field, 
presumed to be the classical depinning field E c i, shows 
the most dramatic increase in CDW current. The lower 
threshold field is interpreted as the Coulomb blockade 
field Ets for soliton nucleation. Figure Qb) shows plots 
of u vs. </>, illustrating the soliton nucleation (0 > tt) 
and classical depinning (0 > C ) instabilities that arise 
as 6 is increased. Figure [4jc) plots u vs. (j) when = tt 
for several values of ue/u p . Figure pld) shows the re- 
sulting phase diagram^ which plots u/tt = cE/(€iEt)) 
vs. ue/u p and allows for variations in e relative to its 
threshold value e\. The diagram illustrates the pinned 
state (6 < tt,ue/u p < 1), a region in which soliton nu- 
cleation occurs (tt < < # c ), and a classical depinning 
region (9 > C ). 

The flat dielectric and other ac response d 13 * 15 * (Fig. [I]) 
and small phase displacements^! below threshold in 
NbSe3 and TaS3 suggest ue/u p « 1 [red arrow in 
Fig. [4] phase diagram] in these samples. Computed (0) 
below threshold^ compares favorably to the reported 
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FIG. 4. (a) Blue bronze I-V curveP 1 in which two threshold 
fields are apparent, (b) Potential energy vs. 0, for increasing 
values of up to the classical depinning instability C . (c) u vs. 
for several ratios ue/u p when = ir. (d) Phase diagram 1 
showing pinned, soliton nucleation, and classically depinned 
states. Red arrow (ue/u p « 1) crosses from the pinned state 
into the soliton nucleation region. The dotted green arrow de- 
picts a system exhibiting both thresholds. Since ue oc 1/e the 
path curves to the left (right) if e increases (decreases) with 
field. Dashed black arrows: classical depinning dominates, as 
suggested by the 4.2-K blue bronze data. 

2° valutP for NbSe 3 provided u E /u p ~ 0.015P Using 
ue/u p = (2ttEt)/E c i , the 48-K blue bronze data-in 
Fig. [4] suggest a similar value of about 0.01. The increase 
in Ets(oc 1/e) with decreasing temperature is readily in- 
terpreted as due to a reduction in e as the normal carrier 
concentration decreases. At 4 K, the normal carriers are 
largely frozen out, resulting in a relatively low e and suf- 
ficiently high Ue/u p for classical depinning to dominate 
[thick black arrows in Fig. [4]. The following sections dis- 
cuss CDW dynamics above threshold and th e issu e of 
quantum coherence, as revealed by CDW ring^SEH and 
other experiments. 



III. TIME-CORRELATED SOLITON 
TUNNELING MODEL 

A basic premise of this paper is that much of the dy- 
namical behavior of CDWs seen in the highest quality 
crystals of NbSe3 and related materials can be under- 
stood by extending the simple picture discussed above. 
These phenomena include narrow-band noise with a fun- 
damental frequency that scales with CDW current and a 
rich spectrum of harmonics, and complete mode- locking 
with an external ac source at high drift frequencies 
(even when much higher than the dielectric relaxation 
frequency, in contradiction with classical predictions^). 
The key to successfully applying such a simple model is to 
accept quantum principles, one of which is Gell-Man's to- 



talitarian principle: 45 "Everything not forbidden is com- 
pulsory." Applied to CDWs the implication is: // CDW 
electrons can tunnel then they must tunnel. Experiments 
to date suggest that CDW condensates behave as sticky 
quantum fluids or deformable quantum solids with dis- 
locations^ 6 rather than massive classical deformable ob- 
jects. 

Hypotheses addressed in this paper include: 1) low 
energy phase soliton dislocations of charge ±2e (or, in 
our view less likely, amplitude solitons of charge ±e/=^"^ 
nucleate above a Coulomb blockade threshold and form 
droplets resembling fluidic domain walls (soliton liquids) , 
where interchain interactions or Josephson coupling be- 
tween chains^ prevent rampant thermal excitations 
2) in the highest quality crystals the nucleation pro- 
cess is best described as coherent Josephson-like tun- 
neling using a modified tunneling matrix element 1 that 
reflects the Zener probability; 3) in these same materi- 
als, the time-evolution of complex order parameters, re- 
sembling probability amplitudes, can be described using 
the Schrodinger equation as an emergent classical equa- 
tion UJ3IJ anc [ 4^ both static (e.g. in ring experiments 
with magnetic flux) and dynamic (ac response) vector 
and scalar potentials can couple to and/or modulate the 
phases of the complex order parameters. 

CDWs are often highly anisotropic, where the dielec- 
tric response, e\\ xi along the chain direction is much 
greater than those, e± yjZ , in the perpendicular direc- 
tions. The degree of anisotropy affects the internal field 
E* generated by a dislocation pair (Fig. [5ja), (b)) and, 
thus, the Coulomb blockade threshold field: E T = E* /2. 
One method of modeling this behavior (using COM- 
SOlJ^) is to rescale the variables along the x-, y-, and z- 
directions by dividing by the relative dielectric constants: 
x' = x/e\\ , y' = x/ej_ y , and z' = z/ej_ z . This is seen start- 
ing with the Maxwell equation: V • D = where (using 
the summation convention): Di = e^e^Ey Here e^- is the 
relative dielectric tensor, which is diagonal with elements 
ey e xx , e ±y e x x, and e ±z e zz , if the axes i,j x,y, 
and z are along the principal crystallographic directions. 
Fig. |5|b) illustrates the rescaled COMSOL simulations 
in 2-D, where the dislocation pair in rescaled coordinates 
looks like a parallel plate capacitor that produces an in- 
ternal field E* = 2e/(2eA c h) = (en c ^)/e, where e = eyeo- 
This is within a factor of - of the ideal value, 2en c ^/e, 
for a fully condensed CDW. Figure [5|c) shows the aggre- 
gation of many 2tt dislocations of charge 2e into fluidic 
soliton droplets that move toward the contacts and al- 
low the bubble of lower energy between them (or 'true 
vacuum,' using the quantum field theory terminology) to 
grow. Other factors that can affect E* and Et include 
gate electrodes in CDW field- and current-effect transis- 
tors EMUl as we n as screening by normal carriers. 

The time-correlated soliton tunneling modelP which 
interprets CDW dynamics above threshold, borrows con- 
cepts from th e the ory of time-correlated single elec- 
tron tunneling! 60 * 61 ! The electrostatic energy parabolas 
of Fig. [3| (also Fig.^a)) are similar to the charging en- 



5 



(a) 0) 




FIG. 5. (a) COMSOL simulation of electrostatic potential 
(red = positive, blue = negative) and field lines for an elec- 
tric dipole consisting of dislocations represented as + and - 
rectangular charge distributions, (b) COMSOL simulation 
for a similar pair with anisotropic dielectric constants, which 
resembles a parallel plate capacitor in rescaled coordinates, 
(c) Aggregation of many dislocations into fluidic domain walls 
or droplets of soliton liquids, between which the bubble of 
lower energy or 'true vacuum' grows as they are driven to- 
ward the contacts by the externally applied field. 



ergies of a small-capacitance tunnel junction. According 
to this model, coherent voltage oscillations, narrow-band 
noise, and ac-dc interference effects come from these 
piecewise parabolic charging energy curves, and not from 
the shape of the periodic pinning potential. The large 
normal carrier concentration in NbSe3 due to incomplete 
Peierls gap formation leads to significant screening by 
normal carriers, which enhances the spatial uniformity of 
the CDW's dielectric response. This explains why such 
highly coherent voltage oscillations, narrow- band no ise 
peaks, and mode locking are often observe d 44 * 62 * 63 * in 
NbSe3 crystals, even though the pinning comes from ran- 
domly distributed impurities EMMl Moreover, the piece- 
wise parabolic curves also explain why the narrow-band 
noise spectra show such a rich array of harmonics. 

The simplest model 1 treats the CDW as a capaci- 
tive Coulomb blockade tunnel junction in parallel with a 
shunt resistor R (Fig. |6jc)) due to normal uncondensed 
electrons. To model dynamics, the 'vacuum angle' is 
related to displacement charge Q between contacts as: 
= 2ttQ/Qo, where Qo = 2eN and N is the number of 
parallel CDW chains. Advancing the phases of all chains 
by 2irn creates multiple pairs of fluidic soliton domain 
walls that quickly reach the contacts. Similar to a capac- 
itive tunnel junction the voltage is: V = (Q — Qq)/2C = 
(Qo/27rC)[6 - 27ra], where C = eA/l. More generally: 
V = (Q /27tC)[0 - (0)], if (0) ^ 27m. The total current 
is: I = I n + I cdw , where I n = (Q /2ttRC)[0 - (</>)] is the 
normal current and I c & w = dQ/dt = (Qo /2tt) d6 jdt is the 
CDW current. (The latter includes capacitive displace- 
ment current but is identical to (Qo/27r)d ((f)) jdt when 
time- aver aged.) Defining uj = 2ttI/Qq and r = RC yields 




FIG. 6. (a) Potential energy vs. for (f> ~ 27m. (b) u vs. (f> 
when — 2tiE/E* > tv(E > Et) as the phases <j>k{x) tunnel 
coherently into the next well via the tunneling matrix element 
T. (c) Time-correlated soliton tunneling model, consisting 
of a normal shunt resistance R in parallel with the CDW, 
represented as a capacitive Coulomb blockade tunnel junction. 



the following equation for the time evolution of 0: 



dO 
~dt 



:U--[0-(4>)\- 



(2) 



Since ((f)) advances in a jerky fashion, Eq. ^ contains 
the elements needed to explain the observed voltage oscil- 
lations, narrow-band noise, etc. Within a unified frame- 
work it allows for at least three mechanisms by which (<p) 
can evolve: (a) coherent Josephson-like tunneling via a 
matrix element T, (b) incoherent tunneling or thermal 
activation of solitons, and (c) classical depinning. Feyn- 
man 31 (vol. Ill, Ch. 21) provides a derivation of coherent 
Josephson tunneling, where the Schrodinger equation is 
viewed as a 'classical' equation to treat wavefunct ion- like 
order parameters coupled by a tunneling matrix element. 
We have developed^ a similar method for the CDW to 
compute (4>(t)) via the coherent tunneling mechanism 
(a). It employs the Schrodinger equation: 



™—^r- = f7o,i^o,i 
at 



1,0 



(3) 



to compute the original and emerging probability ampli- 
tudes ^o(t) & ^i(^) f° r the system to be on branches 
and 1 in Fig. |6|a) (more generally i/jj & V^'+i) when 
coupled by the matrix element T. The model treats the 
amplitudes as complex order parameters: 



^1,0 



where Sq,i = Nq^/N is the fraction of parallel chains on 
the respective branch. Advancing the CDW phases (j)^ (x) 
of many chains by 27r (from one branch to the next in 
Fig. [6| creates lower energy bubbles bounded by droplets 
of microscopic 27r solitons and antisolitons (somewhat de- 
localized as quantum solitons 6 - 4 ^ which form the new flu- 
idic macrostate ipi. 
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The microscopic quantum soliton energy per electron 
pair, A<p, can be estimated from the measured Zener field, 
Eq ~ (A^/^^o e ) 7 typically ~ 10 V/m. Using a phason 
velocity, i/o ~ 3 x 10 3 m/s, yields ~ 5 /ieV, an ex- 
tremely small value. However, the coupled macrostates 
have substantial condensation energies d ue to the many 
(> 10 9 ) interacting parallel CDW chains P^ESl The con- 
densed solitons in the emerging macrostate are thus ef- 
fectively trapped in soliton liquids, preventing thermal 
excitations except across the much larger Peierls gap. 
An analogy is provided by Josephson coupling between 
superimposed macrostates in 2-band superconductors,^ 
where thermal excitations only occur across either BCS 
energy gap regardless of the energy difference between 
macrostates. One can also view bubbles of the CDW 
chains escaping out of the metastable well (Fig. [6jb)) 
as being analogous to superfluid helium atoms quantum 
mechanically creeping out of a container. If the con- 
tainer rim is, for example, d ~ 1 cm above the liquid sur- 
face, then the gravitational barrier per atom is mgd ~ 
4 neV, which is small compared to kT even at 1 mK. 
Nevertheless, the helium atoms remain trapped in the 
superfluid, prevented by the condensation energy from 
thermally hopping out of the container even though they 
quantum mechanically creep over the rim in a collective 
fashion. 

The driving force F is the energy difference per unit 
length after one branch crosses another in Fig. j6^a)P^ 
Using t he an alogy to pair production^ and following 
Bardeen j 66 ! 67 ' the tunneling matrix element T is esti- 
mated to be: 



(a) 



2.0 



(b) 



T(F) = -4FAexp [-F /F] , 



(5) 



where Fq (A^)/(7ico) and A is defined in ref. 1 . Fig. [7|a) 
compares the simulations with measured voltage oscil- 
lations^ of NbSe3 for rectangular current pulses. Ex- 
cept for the increasing pulse amplitudes, the same pa- 
rameters are used for the entire family of theoretical 
plots (solid lines), which show unprecedented quantita- 
tive agreement with experiment. The model correctly 
captures the progression of non-sinusoidal shapes, rang- 
ing from rounded backward sawtooth behavior for the 
9.90-/iA current pulse to more symmetrical oscillations 
for higher pulse amplitudes. The inset to Fig.[7ja) shows 
the CDW current (I c dw — I ~ In) vs - time correspond- 
ing to the 10.89-/iA pulse. This plot: 1) shows that 
a large fraction of the CDW current is oscillatory, and 
2) captures the 'flowing,' rather than abrupt tunneling, 
aspect of quantum transport. The I-V and differen- 
tial resistance curves are computed 1 by averaging over 
several cycles, with results shown in Figs. [T^b)-(d). A 
range of behaviors are captured, ranging from rounded 
Zener-like behavior to more linear I-V curves and dV/dl 
curves with negative dips or wings, as seen in NbSe3 crys- 
tals with fewer impurities. 44 The theoretical plots show 
outstanding quantitative agreement with experiment in 
Fig. |7|;d). 
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FIG. 7. (a) Theoretical (solid lines) vs. experimental (dotted 
lines^) voltage oscillations (bottom to top, offset by 0, 0.25, 
0.5, and 0.75 mV) of an NbSe3 crystal at 52 K for current 
pulse amplitudes: 9.90 fiA (black), 10.89 fiA (red), 11.49 fiA 
(green), & 11.88 fiA (blue). Inset. CDW current, I - In, 
vs. time for the 10.89 fiA pulse, (b) Simulated CDW current 
vs. field for several go = Fo/2eE*. Dotted lines: Bardeen's 
modified Zener function.^ (c) Simulated R = dV/dl vs. cur- 
rent for several go, where R n is the normal resistance below 
threshold, (d) Theoretical (solid lines) vs. experimental (dot- 
ted lines) dV/dl vs. current for NbSe3 (see^for parameters). 



IV. COUPLING OF ORDER PARAMETER 
PHASES TO VECTOR AND SCALAR 
POTENTIALS 

The simulations of dc transport and rectangular cur- 
rent pulses, discussed above, fix the phase difference 
S = 5j+i — Sj at 7r/2 in Eq. (p])P which yields the max- 
imum current in the Josephson current-phase relation. 
This section discusses coupling of static (magnetic field) 
and dynamic (ac electric field) vector and scalar po- 
tentials to the phases Sj in order to: a) inter pret t he 
h/2e quantum interference effects in CDW rings J^^b) 
better understand mixing and other ac response experi- 
ment EMII] previously interpreted using photon- assisted 
tunneling (PAT) theory}^ and c) interpret new large am- 
plitude ac experiments presented here. The ability to 
vary both magnitudes and phases of the macrostate am- 
plitudes could eventually set the stage for development of 
future quantum computing devices, while better under- 
standing of the ring experiments could enable new types 
of magnetic sensors. 

In the CDW ring experiments^^ a static magnetic 
vector potential couples to the phases Sj and leads to 
quantum interference between the amplitudes traversing 
the two branches of the ring (nucleated quantum solitons 
interfering with themselves). This can be visualized in 
terms of an extra phase shift x affecting the tunneling 
matrix elements of a two-domain model: 
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T a , b ^Texp[±i X /2], (6) 
one domain for each path, a or b, along the ring. Here: 



A-dr = 2ir [*/$ ] ■ 



(7) 



where $o = h/q and # is either e or 2e. Summing 
the amplitudes then yields a modulation proportional to 
\2T cos[$/$o] | • This simple two-domain picture gives the 
correct period h/2e for the A-B oscillations provided we 
take q = 2e. However, it is an oversimplification com- 
pared to the reported ~ 10% modulation and rather dis- 
ordered be havio r in the A-B oscillations in the actual 
CDW ring d 1Q | 11 J , suggesting the need to include many 
CDW domains with some degree of disorder. More- 
over, ac response experiments, discussed below, suggest 
a rather short tunneling distance. This further indicates 
the need to incorporate multiple domains, which can be 
modeled as a network of many tunnel junctions in series. 

The phases Si(t) of the macrostate order parameters in 
Eq. Q can be modulated by an oscillatory field that tem- 
porally evolves the scalar and/or vector potentials. The 
theory of photon- assisted tunneling (PAT) 69 enables pre- 
dictions of tunnel junction response to combined dc and 
ac signals based on its dc current- volt age (I-V) charac- 
teristic. Oscillatory voltages modulate the relative en- 
ergies and phases of wavefunctions on opposite sides of 
the tunnel junction. This generates various combinations 
of Bessel functions in the predicted responses, which re- 
duce to finite differences of the I-V curves in the small- 
signal limit. A modification of PAT theory was previ- 
ously adapted to interpret mixing and other CDW ac 
response experiments^"^ on TaS3 and NbSe3. These 
experiments show good agreement with PAT theory for 
small-amplitude signals! 13 " 15 ! Mixing experiments apply 
a signal of the form: 



V(t) = Vdc + V\ COSUJit + Vi COSCc^ 



(8) 



and measure an induced response (e.g. with a lock-in am- 
plifier): SI(t) = SIo cos[<jJot+(p]. The difference frequency 
is: ujq = |co>2 — oo i\ for direct mixing and uoq = \oo2 — 2ooi\ 
for harmonic mixing. At low frequencies and amplitudes, 
the harmonic mixing response vs. bias voltage Vd c is pro- 
portional to the third derivative of the dc I- V curve: 



SIo(V dc ) = -V?V 2 



d 3 I 



dc 



dV 3 



(9) 



v=v dc 



At finite frequencies, the third derivative gets replaced 
by a third finite difference^ with a step size propor- 
tional to frequency but has a similar, albeit broadened, 
bias dependence. The harmonic mixing response at zero 
dc bias voltage becomes significant for frequencies ooo/2tt 



of about 1 MHz and greater, and is found to be bias- 
independent below threshold.^ 

Rather different behavior emerges when the dc bias 
voltage Vdc hi Eq. ([8| is replaced by a large amplitude ac 
'bias' voltage, 



V dc -> V dc cos UJt, 



(10) 



and the harmonic mixing response SIo is plotted vs. V ac . 
When uo is small, since harmonic mixing is an even func- 
tion of dc bias, SIo(V ac ) is just the time-averaged re- 
sponse vs. bias voltage (Wo (Vbias (£)))> which resem- 
bles a washed out third derivative. When uo/2tt reaches 
about 50 kHz or higher, however, SIo(V ac ) resembles 
SIo (Vdc — Vt) with an apparent threshold voltage col- 
lapsed to the origin 15 , possibly because of capacitive cou- 
pling due to the extremely high CDW dielectric response. 

The most interesting behavior is expected to occur 
when uo is at megahertz frequencies and higher. In the 
absence of coupling T, a state of modulated energy Eq,i 
corresponding to branch or 1 in Fig. [6] would evolve as: 



^o,i (t) = ^o,i (0) exp 



/ft) f dt'Eo,i(t') 
Jo 



(ii) 



Taking the charge to be 2e, the voltage Vi(t) = Vg cosoot 
across a small domain of length i modulates the energy 
E\(i) of state 1 relative to Eq, as: AE(t) = 2eVecosoot. 
State ipi then evolves relative to as: 

ipi (t) = ipi (0) exp [— iz sin oot] 

oo 

= ^o(0) ^2 J n(z) exp [-inuot] , (12) 



n= — oo 



where J n (z) are Bessel functions and z = (2eVi)/hoo. 
This effectively splits up the ipi amplitude into many, 



tfn = Jn(z)lpl, 



(13) 



of virtual energy E n = nftoo. 'Turning on' the tunnel- 
ing matrix element T enables it to couple states ifjo & 
%jj' n of equal energy, any negative energy difference being 
balanced by the soliton pair energy. 

Equation (13) thus captures the essence of photon- 



assisted tunneling, where an initially occupied state can 
tunnel into an unoccupied state of equal virtual en- 
ergy. Recalling the relation between harmonic mixing 
and dc bias voltage (E q. j9| and finite difference forms^), 
following PAT theoryp^ESJ and noting that J- n (x) = 
( — l) n J n (x), the harmonic mixing response vs. total volt- 
age amplitude V ac between contacts would then be ex- 
pected to be given by: 



SI (V ac ) = — ) SI (Vdc = 0) 
* aoo 1 



n—1 x 7 



naoo) (14) 
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where Vb = Vdc — Vt due to the collapsed effective I-V 
curve^ at finite frequencies. The amplitudes and fre- 
quencies, Vi, V2 , cji, and cj2> of the signals inducing the 
harmonic mixing response are fixed in these experiments. 



The scaling parameter a in Eq. (14) depends on the 
distance L between contacts and an effective scaling 
length which relates the energy acquired by a parti- 
cle of charge e* in an electric field E to a quantum of 
energy Tiuj: e*Vi = e*E£ o Tiuj. Previously,^ the effec- 
tive charge was assumed to be a reduced by the Frohlich 
mass ratio, e* ~ [tu/Mf] e ~ 10 _3 e, which yielded values 
of £ in the range 1.6-22 ^ m LLU M otivated by the recent 
CDW A-B ring experiments,^^ here we take the effec- 
tive charge to be: e* = 2e, which reduces the estimated 
values of £ into the nanometer range. In general, the 
nature of Zener-like tunneling through a tilted soliton 
gap may yield some degree of frequency- and/or field- 
dependence of £. Using V ac = (L/£)Vi, one obtains the 
following scaling parameter: 



L h 



(15) 



Due to the properties of Bessel functions, the Jq (V ac /auj) 
term in Eq. ( [14] ) should initially dominate for small am- 
plitudes, Vac, while the remaining terms may become sig- 
nificant for larger V ac . Defining SI m = SIo(V ac = 0), re- 
stricting A n (uo) = 25Io(nau)/5I m to be real, and keeping 
a finite number, iV, of terms, the normalized theoretical 
harmonic mixing response can be approximated as: 



SIo_ 

Sim, 



T 2 



Vac 



N 

£ 

n=l 



Vac 

auj 



(16) 



Figures [8] and [9] show plots of normalized harmonic 
mixing responses 1 52o I / 1 I vs - ac bi as amplitude V ac 
for single crystals of TaS3 and NbSe3, as compared to 
Eq. (16). Figure [8ja) shows measured harmonic mixing 
responses of a TaS3 crystal (L = 0.1 mm) with a 5 mV 
threshold voltage at 180 K, for three different ac bias 
frequencies uj. Figure j8^b) shows theoretical plots using 
Eq. (16) and the parameters in Table [ij For this sample 
the effective scaling distance £, estimated using Eq. ( 15 ) 
from the parameter a, is found to be in the range 8- 
15 nm, or several CDW wavelengths. The extremely 
small soliton energy gap per electron pair enables this 
distance to be longer than one would normally encounter 
in an ordinary tunnel junction. 



TABLE I. Eq. |l6} parameters used for the Fig. gb) theo- 
retical plots. 



cj/2tt (MHz) 


a(V-s) 




A 1 


A 2 


A 3 


A A 


40 


5.5 x 10" 


12 


1.75 


0.25 


0.42 


-1.00 


80 


4.2 x 10" 


12 


1.80 


1.40 


0.50 


-0.60 


160 


3.2 x 10" 


12 


2.00 


0.70 


1.57 


1.44 
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V ac (mV) 



10 20 30 40 50 60 
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FIG. 8. (a) Normalized magnitude of harmonic mixing re- 
sponse ( cji/2tt = 5 MHz, cj 2 /2tt = 14 MHz, cj /2tt = 4 MHz) 
of a 0.1-mm long TaS3 crystal, with a dc threshold Vt — 
5 mV, vs. ac bias amplitude V ac at three different frequencies 
(jj/2tt at 180 K. (b) Theoretical plots using Eq. (16) and the 
parameters shown in Table [I] 



Figure [9] shows measured harmonic mixing responses 
of an NbSe3 crystal (L = 5 mm) at 120 K, for several 
ac bias frequencies uj. Figure [9Fa) shows experimen- 
tal plots, while Fig. |9jb) shows theoretical plots using 
Eq. JT6] ) and the parameters in Table |n[ Figures [9jc) 
andpMd) directly compare experiment with theory for ac 
bias frequencies of 4 MHz and 8 MHz. For this sample 
the effective scaling length estimated using Eq. (15) 



from the scaling parameter a, is found to be 1.5 nm, or 
slightly greater than one CDW wavelength. Here, the 
shorter length £ may reflect a reduced effective mean free 
path length for the quantum solitons due the incomplete 
Peierls gap and large number of uncondensed normal car- 
riers in NbSe3. It will be interesting, in future studies, 
to determine whether the CDW wavelength represents a 
lower bound on £. 



TABLE II. Eq. ( 16 ) parameters used for the Fig. [^theoretical 
plots. 



a = 1.1 x 10~ 9 (V-s) 


Ai 


A 2 


A 3 


A 4 


4 MHz 


0.40 


0.65 


0.15 


0.55 


5 MHz 


0.45 


0.70 


-0.30 


0.10 


6 MHz 


0.30 


1.00 


-1.40 


1.90 


7 MHz 


0.20 


1.10 


-2.00 


3.50 


8 MHz 


0.20 


1.40 


-1.65 


0.70 


10 MHz 


0.15 


1.90 


-1.90 


-1.95 


15 MHz 


0.01 


3.00 


3.00 


1.00 



The experiments reported here, as well as earlier mix- 
ing experiments at temperatures sometimes exceeding 
200 Kp^HH] are consistent with the idea that oscilla- 
tory electric potentials modulate the phases of classically 
robust order parameters resembling macroscopic wave- 
functions. Moreover, the experimental results are consis- 
tent with those of the CDW ring experiments j 10 ! 11 ! which 
demonstrate a significant degree of CDW quantum co- 
herence. Collectively, the experiments support the hy- 
pothesis that either a vector or scalar potential couples 
to order parameter phases of CDW soliton condensates, 
and in some cases can lead to quantum interference. Fur- 
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FIG. 9. (a) Normalized magnitude of harmonic mixing re- 
sponse ( cji/2tt = 1 MHz, (jJ 2 /2tv = 2.8 MHz, oj /2ty = 
800 kHz) of an NbSe3 crystal (L = 5 mm), vs. ac bias ampli- 
tude Vac at seve ral frequencies cj/27t at 120 K. (b) Theoretical 
plots using Eq. ( 16 ) and the parameters shown in Table [Tl] (c) 
and (d) Direct comparisons between theory and experiment 
for frequencies oj/2tt of 4 MHz and 8 MHz. 



ther experimental and theoretical studies are warranted 
to enable the eventual development of a microscopic de- 
scription of CDW transport. 



V. DISCUSSION AND CONCLUSION 

CDW transport is one of the few known cases of cor- 
related transport of macroscopic numbers of electrons - 
the only known example of large-scale collective electron 
transport at human body temperaturesP This paper is 
highly transformative in that it challenges the classical 
sliding CDW paradigm that has dominated the field for 
over thirty years. Nevertheless, the quantum ideas dis- 
cussed here can hardly be regarded as speculative. The 
evidence supporting quantum theory is so enormous, it 
can be considered a proven fact that electrons and all 
other known particles behave quantum mechanically. In 
2000, the 100-year anniversary of Planck's black-body ra- 
diation paper, 72 Kleppner and J ackiw 7 -^ pointed out that: 
"Quantum theory is the most precisely tested and most 
successful theory in the history of science." Since then, 
aspects of quantum theory (the Pauli principle 7 -^) have 
been confirmed to within an accuracy of 6 x 10 -29 . 

The classical behavior one observes on the macroscopic 



scale depends on the system and emerges from the be- 
havior of large numbers of entangled quantum particles 
exhibiting wave-particle duality. The Schrodinger equa- 
tion can be regarded as the 'classical' equation for super- 
conducting condensates coupled through a thin insulator 
by Josephson tunneling (vol. Ill, Ch. 21 oP^). Similarly, 
the time-correlated soliton tunneling model discussed 
here treats the Schrodinger equation as an emergent clas- 
sical equation describing Josephson-coupled fluidic CDW 
macrostates. The simulations yield unprecedented quan- 
titative agreement with coherent voltage oscillations and 
I- V characteristics of NbSe3 and also provide a natural 
interpretation for the quantum interference seen in the 
CDW ring experiment J 10 * 11 * and more complex interest- 
ing behavior seen in CDW harmonic mixing response. 

Any further progress in understanding of CDW trans- 
port will require the scientific community to accept the 
fact that the CDW electron-phonon condensate behaves 
according to laws of quantum physics - the same quan- 
tum principles that govern every other system of parti- 
cles in the universe. It is not necessarily true, a priori, 
that quantum principles are consistent with the current 
dogma - that CDW electrons classically "slide" according 
to Aristotle's linear velocity-force relation. Addressing 
the quantum behavior of CDWs, perhaps culminating in 
a microscopic theory of CDW transport and dynamics, 
would have enormous impact on this important branch 
of condensed matter physics. Additional areas of broad 
impact potentially include the boundary between CDWs 
and superconductivity, correlated electron-ion transport 
in biological systems, tunneling and 'false vacuum decay' 
in quantum cosmology, a formally similar — it instabil- 
ity for spontaneous CP violation,^ and a deeper under- 
standing of quantum theory and the emergence of classi- 
cal reality at the macroscopic scale. Observation of quan- 
tum effects in NbS3, which undergoes a Peierls transition 
well above room temperature,^ would potentially lead to 
new devices such as magnetic sensors operating at room 
temperature. Understanding of the quantum behavior of 
solitons could lead to topologically robust (against deco- 
herence) forms of quantum information processing, which 
would have major technological significance. 
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